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We investigate, within a local density functional theory formalism, the interactions between like- 
charged polyions immersed in a confined electrolyte. We obtain a simple condition for a repulsive 
f*^ , effective pair potential, that can be related to the thermodynamic stability criterion of the uncharged 

counterpart of microscopic species constituting the electrolyte. Under the same condition, the 
' phenomenon of charge inversion (over-charging), where the polyion bare charge is over-screened by 

its electric double layer, is shown to be impossible. These results hold beyond standard mean-field 
theories (such as Poisson-Boltzmann or Modified Poisson-Boltzmann approaches). 
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PACS numbers: 05.70.Np, 64.10+h, 82.60Lf, 82.70.Dd 



Electrostatic forces play a key role in determining the stability, phase and structural properties of colloids or 
macroions suspended in polar solvents (usually aqueous media). Such suspensions are ubiquitous in many biologi- 
cal and technologically important systems, from superabsorbants or water-soluble paints to DNA solutions. Their 
behavior is however far from being well understood, despite an intense theoretical effort over the last fifty years. 

Upon dispersion in water, mesoscopic colloidal polyions bearing ionizable groups gain a (large) bare charge, and 
strongly attract counterions while repelling coions, thus giving rise to electric double layers characterized by strong 
inhomogencitics in the local density of these microions in the vicinity of the polyions. From the electrostatic part of the 
traditional Derjaguin-Landau-Verwey-Overbeek (DLVO) theory H, the effective interactions (solvent and microion 
mediated) between two like-charged colloidal particles are expected to be repulsive whereas experiments || and 
numerical simulations Q provide evidence for attraction, particularly within confined geometries (e.g. when the two 
O \ macroions are close to a charged wall or between two glass plates). 

We consider a mixture of N microions, where species a has charge number z a and local density c Q (r). Disregarding 
\ the molecular nature of the solvent, considered to be a confined electrolyte of dielectric permittivity e, the effective 
^ ■ interactions between two like-charge colloids are then analyzed within a generic local density functional theory. We 
0^ ' show that the convexity of the free energy density with respect to the microions densities is a sufficient condition for 
having repulsive effective interactions and for preventing charge inversion. The physical implications of this result 
are considered in the final discussion (in particular the repulsion "often" follows from the thermodynamic stability 
criterion of the corresponding uncharged mixture of micro-species). The proof makes use of the results established 
in Refs. [^[]7); and considers a geometry which encompasses many cases of experimental relevance: the confining 
region TZ is a cylinder of arbitrary cross section and length 2L (in the limit of large L) while the medium TZ 1 outside 
TZ is a dielectric continuum of permittivity e' with the possibility of a uniform density of surface charges a on the 
boundary dTZ. The standard Neuman and Dirichlet boundary conditions (constant normal electric field and constant 
potential respectively) belong to the above class, and the dielectric medium in TZ' may contain an electrolyte solution 
(this possibility was not considered in Refs We assume global electroneutrality in both TZ and TZ' and the two 

macroions we shall focus on may be of arbitrary shape provided the electrostatic potential possesses mirror symmetry 
with respect to a plane Oxy between the macroions. Our approach holds irrespective of the specific boundary 
\ conditions to be applied on the macroions (constant charge or constant potential) and is independent of the sign of 
the surface charge a. 

Omitting, for the sake of simplicity, their temperature dependence, the local || density free energy functional || 
considered here are of the form 



H{c a }) = f f({c a (r)})dr + \ [ p c (r)G(r,r') Pc (r')drdr', (1) 

where the Helmholtz free energy density / is a function of all TV microions densities and the local charge density 
includes the contributions from the macroions q M (r) as well as from the microions: 

N 

Pc(r) = J2z a e c a (r) + q M (r). (2) 
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In Eq. ([!]), G(r, r') denotes the Green's function such that the electrostatic potential ^>(r), which is a solution of 
Poisson's equation V 2 V> = — (47r/e)/9 c (r) with the required boundary conditions on the surface dlZ, can be cast in the 
form 

^(r)= / p c (r')G(r,r')dr'. (3) 



G reduces to the usual Coulomb potential in the absence of confinement. The density functional formulation provided 
by Eq. (Q) encompasses a large class of modelizations, including a) the standard Poisson Boltzmann theory when 
the ideal gas expression / = ^ Q fcTc Q ,[ln(c Q ,A^) — 1] is chosen for the free energy density, b) modified Poisson 
Boltzmann theories derived to account for steric effects or more general non-electrostatic interactions JllJ , and c ) 
recent improvements over standard mean-field approximations incorporating ion-ion correlations in the form of a One 
Component Plasma correction |L2| . This last case involves a free energy density which does not only depend on the 
local densities of microions but also on the elementary charge e, and will be examined in the concluding discussion. In 
general, the present analysis allows to incorporate correlation or non mean-field contributions, provided these effects 
translate into a local free energy density. 

Minimization of the functional (II]) subject to the normalization constraints J^c a (r)dr = N a yields the implicit 
relation between the densities c Q (r) and the electrostatic potential: 

df 

+ e z Q V>(r) = AC l<a<iV (4) 

UCot 

where (tx* is the (electro) chemical potential of species a. The effective force acting on a colloid is obtained from 
integration of the generalized Maxwell stress tensor II over the surface S of the macroion: 

U-ndS, (5) 

s 

where the unit vector n points outward from the surface of integration. The local stress tensor follows from the 
computation of the reversible work required to locally deform the shape of the confining cylinder, with the result 



n 



P + ~ Vi/j ■ ViJj 



eVipfgiVip, (6) 



where the osmotic pressure P({c a }), which depends on position via the local densities c a (r), is the pressure that 
would be found in the medium in the absence of an electric field, with the same densities c a (r) fl3|| , and can be 
written as the Legendre transform of the free energy density: 

N 

^(K}) = -/(K}) + E c «af- (?) 

a— 1 a 

For any stress tensor of the form (||), it has been shown in Refs. that the effective force F given by Eq. (|) has 
a component F z along the axis of the cylinder that can be cast in the form 



Oxy 



dP 



dxdy 

+ | / (V^ =0 - VV^ =i ) 2 dxdy, (8) 

z JOxy 

where e denotes s (respectively e') for the part of the symmetry plane, denoted Oxy, belonging to 1Z (respectively 
TV). In Equation (|J), it is understood that the pressure depends on the potential via Eqs. (0) and (Q), the latter 
giving sense to the notation dP/dip. The subsequent analysis is devoted to the proof that P is a convex- up function 
of ij} provided that the free energy density is itself convex- up with respect to the densities {c a }. Remembering the 
orientation convention chosen in Refs. this in turn establishes the repulsive nature of the effective interactions 

(F z > 0). From the direct differentiation of Eq. (|?|) we get 
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where use was made of the stationary condition (Q) . The relation (Q) can be recovered from the mechanical equilibrium 
condition of a fluid element of microions: the balance between the electric force and the osmotic constraint acting on 
such a fluid element located at point r can indeed be written 



VP 

dP 

dip 



p c (r) Vip 



Vip = 0. 



(10) 

(11) 



Equation (||) equivalently implies that the tensor II is divergence free, as already invoked in Ref. to obtain the 
osmotic equation of state in the specific case of a Modified Poisson Boltzmann theory. 
The second derivative of the osmotic pressure, 



d 2 P 



dc a 



dip 2 ^ " dtp 

a 

involves the quantities dc a /dip which may be determined from the condition (^): 



(12) 



JV 



\ - d 2 f dc p 

^ dc a dcp dip 



-ez a 



(13) 



Assuming the free energy density to be a convex- up function of the {c a } implies that the matrix H with elements 

d 2 f 



H, 



a/3 



dc a dc/3 

is positive definite. Its inverse if -1 is consequently positive definite from which we can invert relation (p 



dc a 
dip 



N 
(8=1 



>a/3 



Zp. 



Upon substitution of Eq. ( |l5| ) into (|T^) we finally obtain 

d 2 P 



dip 2 



= e 2 E( H_1 )a/3^^ >0 



(14) 



(15) 



(16) 



a, 



The pressure P is therefore a convex-up function of the electrostatic potential, which establishes the aforementioned 
connection between the convexity of / and the sign of the polyion-polyion effective interactions. Note that the non 
convexity of the free energy density is a necessary though insufficient condition for the existence of attraction within 
a local density approximation formalism. 

We now address the question of charge reversal, where the nominal charge of a bare polyion can be over-screened by 
its condensed counterions. This interesting phenomenon (occurring for instance in DNA salt systems jli]]) requires, 
from Gauss's theorem, the vanishing of the normal component of the electric field (— n • Vip) on a closed surface with 
normal n around the polyion. Taking the gradient of Eq. (|j) then implies, without any symmetry assumption 



N 



d 2 f 

deadep 



[n • Vc p ] = 0, 



(17) 



hence the existence of a vanishing eigenvalue [with eigenvector (n • Vc ( g)/3 = i...Ar] of the stability matrix H, which is 
impossible under the assumption of positive definiteness. 

Within the framework of Poisson Boltzmann theory, the repulsive nature of the effective electrostatic force between 
two similar colloidal bodies (^0] and the impossibility of charge reversal immediately follows from the previous 
analysis. However, if for the sake of analytical tractability and in the spirit of Ref. the ideal free energy density 
is further Taylor expanded up to order Af around the mean densities {c Q }, we obtain 



x ^~ 2 

° a i=0 



M-l 



Sg(3 



(18) 
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The truncated Hessian H is therefore positive definite and gives a repulsive force for even values of Af (as in the 
non-linear Poisson-Boltzmann theory), while it can be non-positive definite and the force attractive for odd orders Af 
provided the local density of one counterion species exceeds twice its mean density. This explains the repulsion within 
linearized Poisson-Boltzmann theory [[l6f (where linearization amounts to taking Af = 2 ||), and the attraction seen 
in Ref. Jl5| with an expansion pushed one order further and truncated after third order to account for the triplet 
interactions between two polyions and a charged wall. 

We have shown that effective attractive pair potentials and over-charging are both ruled out under the assumption 
of positive definiteness for the stability matrix H. This is a necessary requirement when the free energy density does 
not depend on the elementary charge e. The key ingredient here is the existence of a thermodynamically stable neutral 
mixture of micro-species described by the same density f{c a } as its charged counterpart. On the other hand, the 
inclusion of correlations and/or fluctuations corrections is believed to be essential in order to describe the appearance 
of attractive interactions. If their correction to the mean-field electrostatic term J p c Gp c in the right hand side of Eq. 
(|l|) is taken into account, and translated into an e-dependent free energy density (T^JT^], the thermodynamic stability 
condition is now the positive definiteness of the integral operator whose kernel is defined by 

' hJ 6(t-t>) +e 2 z a z G(v,v'), (19) 



dc a dcp 

and it is in general impossible to find a locally neutral mixture described by the same density /. H need not therefore 
be positive definite. However, non convex-up densities / such as the correction put forward by Stevens and Robbins 
fl7f by extrapolation of the Monte Carlo data for the homogeneous One Component Plasma lead to an instability at 
high densities (the so called structuring catastrophe), and cannot be used in a thermodynamically stable way within 
a local density functional theory [fl2||T§| l . The recent attempt of Barbosa et al. Jl2| to circumvent this instability gives 
rise to a convex-up /, with results reproducing satisfactorily the ionic correlations present in Molecular Dynamics 
simulations. The present work indicates that the effective interactions are then necessarily repulsive and charge 
inversion impossible, regardless of the analytical complexity of the formalism. 

To conclude, repulsive effective pair potentials seem generic within local density approximations. The description of 
attraction and over-charging thus requires involved density functional theories, and it is not sufficient to incorporate 
correlation effects in a local formulation. Our analysis points to the importance of non-local effects, that can be 
accounted for by weighted density functionals ]T^Jlq| , or simpler approaches relying on non-electrostatic depletion 
interactions JT9| . 
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and D. Rowan for helpful discussions. 
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